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Dynamical horizons are considered in full, non-linear general relativity. Expressions of fluxes of 
energy and angular momentum carried by gravitational waves across these horizons are obtained. 
Fluxes are local, the energy flux is positive and change in the horizon area is related to these fluxes. 
The flux formulae also give rise to balance laws analogous to the ones obtained by Bondi and Sachs 
at null infinity and provide generalizations of the first and second laws of black hole mechanics. 
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Black holes which are themselves in equilibrium but in 
possibly time-dependent space-times can be modelled by 
isolated horizons [1]. Over the past three years, proper- 
ties of isolated horizons were studied in detail. In par- 
ticular, the framework enabled one to extend the laws of 
black hole mechanics [2, 3] and has been used to extract 
physics from initial data of widely separated black holes 
[4] and from numerical simulations of the final phases of 
black hole collisions [1, 5]. The purpose of this letter is 
to outline a generalization of these ideas to fully dynami- 
cal situations in which matter and gravitational radiation 
can fall into black holes. 

Our definition of a dynamical horizon is geared to prac- 
tical applications in astrophysical situations, particularly 
those considered in numerical relativity. 



I. DEFINITION AND NOTATION 

Definition: A smooth, three-dimensional, space-like 
sub-manifold H in a space-time is said to be a dynamical 
horizon if it is foliated by a preferred family of 2-spheres 
such that, on each leaf S, the expansion 9m of a null 
normal £ a vanishes and the expansion 9i n ) of the other 
null normal n a is strictly negative. 

Thus, a dynamical horizon H is a 3- manifold which 
is foliated by marginally trapped 2-spheres. As shown 
below, the area of these 2-spheres necessarily increases. 
An example is provided by continuous segments of world 
tubes of apparent horizons that naturally arise in numer- 
ical evolutions of black holes. In contrast to event hori- 
zons, dynamical horizons can be located quasi-locally; 
knowledge of the full space-time is not required. The 
condition that H be space-like is implied by a stronger 
but physically reasonable restriction that the derivative 
of 9(t) along n a be negative [6]. Finally, the requirement 
that the leaves be topologically S 2 can be replaced by 
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the weaker condition that they be compact. One can 
show that the topology of S is necessarily S 2 if the flux 
of matter or gravitational energy across H is non-zero. If 
these fluxes were to vanish identically, H would become 
isolated and replaced by a null, non-expanding horizon 
[2]- 

Dynamical horizons are closely related to Hayward's 
trapping horizons [6]. However, while the definition of 
trapping horizons imposes a condition on the derivative 
of off H, our conditions refer only to geometric quan- 
tities which are intrinsically defined on H. But in cases of 
physical interest, the additional condition would be satis- 
fied and dynamical horizons will be future, outer trapping 
horizons. Nonetheless, our analysis and results differ con- 
siderably from those of Hayward's. While his framework 
is based on a 2+2 decomposition, ours will be based on 
the ADM 3+1 decomposition. Our discussion includes 
angular momentum, our flux formulae are new and our 
generalization of black hole mechanics is different. Our 
analysis is geared to providing tools to extract physics 
and perform checks on numerical simulations of dynam- 
ical black holes. Therefore we will restrict ourselves to 
dynamical horizons with zero charge. 

Let us begin by fixing notation. Let r a be the unit 
time-like normal to H and denote by V a the space-time 
derivative operator. The metric and extrinsic curvature 
of H are denoted by q a t and K a f, :— q a c qb d ^cTd respec- 
tively; D a is the derivative operator on H compatible 
with q a b and 1Z a b its Ricci tensor. Leaves of the preferred 
foliation of H will be called cross- sections of H. The unit 
space-like vector orthogonal to S and tangent to H is de- 
noted by r a . Quantities intrinsic to S will be generally 
written with a tilde. Thus, the two-metric on S is q a b, the 
extrinsic curvature of S C H is K a b '■= Q a °% DcTd, the 
derivative operator on (S, q a b) is D a and its Ricci tensor 
is IZab- Finally, we will fix the rescaling freedom in the 
choice of null normals via l a := r a +r a andn a := f a — r a . 

We first note an immediate consequence of the defini- 
tion. Since 6m = and 9i n ) < 0, it follows that K > 0. 
Hence the area as of S increases monotonically along r a . 
Thus the second law of black hole mechanics holds on H . 
We will obtain an explicit expression for the change of 
area in part III. 
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Our main analysis is based on the fact that, since H 
is a space-like surface, the Cauchy data {q a b,K ab ) on H 
must satisfy the usual scalar and vector constraints 

H s := K + K 2 - K ab K ab = \^GT ab r a r b (1) 
Hy := D b (K ab - Kq ab ) = 8irGT bc ? c q a b . (2) 

We will often focus our attention on a portion AH C H 
which is bounded by two cross-sections Si and S 2 - 



II. ANGULAR MOMENTUM 

The angular momentum analysis is rather straight for- 
ward and is, in fact, applicable to an arbitrary space-like 
hypersurface. Fix any vector field ip a on H which is tan- 
gential to the cross-sections of H. Contract ip a with both 
sides of eqn. (2). Integrate the resulting equation over 
the region AH C H , perform an integration by parts and 
use the identity L v q ab = 2D^ a ip^ to obtain 

^i K ^ a?bd2V -^i K ^ bd2V 

^ Aff (^^li rt ^) rf3y (3) 

where P ab := K ab - Kq ab . It is natural to identify the 
surface integrals with the generalized angular momentum 
associated with those surfaces and set: 

J s ] --^cj s K ^ a r b d 2 V (4) 

where we have chosen the overall sign to ensure compat- 
ibility with conventions normally used in the asymptoti- 
cally flat context. The term 'generalized' emphasizes the 
fact that the vector field ip a need not be an axial Killing 
field even on S; it only has to be tangential to our cross- 
sections. 

The flux of this angular momentum due to matter 
fields and gravitational waves are respectively 

Ji v) = - f T ab T a v »d 3 V, (5) 

JAH 

J™ = -^af^^V, (6) 
and we get the balance equation 

4 v) -j[ v) =jM+jM. (7) 

As expected, if (p a is a Killing vector of the three-metric 
q a b, then the gravitational flux vanishes: = 0. 

It is convenient to introduce the angular momentum 
current j v := —K ab ip a r b so that (4) becomes — 
(8TrG)-i$ s j*d*V. 



III. ENERGY FLUXES AND AREA BALANCE 

As is usual in general relativity, the notion of energy is 
tied to a choice of a vector field. Here, we will consider 
vector fields £ a = N£ a where the lapse N is constructed 
as follows. Let r be a radial coordinate on H defined such 
that the cross-sections of H are level surfaces of r. Then 
r a oc D a r. It turns out that in order to get the balance 
law for energy, we must tie our lapse functions N to radial 
coordinates such that D a r = N r r a . (Since £ a = N£ a , 
as usual the term 'lapse' refers to space-time evolution; 
not to 'evolution' along f a .) Thus each r determines a 
permissible lapse function N r . If we use a different radial 
coordinate r' , then the lapse is rescaled according to the 
relation 

fir' 

N r ,=N r —. (8) 
ar 

Thus, although the lapse itself will in general be a func- 
tion of all three coordinates on H, the relative factor be- 
tween any two permissible lapses can be a function only 
of r. We denote the resulting permissible vector fields by 
£" r ) := N r £ a . Recall that, on an isolated horizon, physi- 
cal fields are time independent and null normals can be 
rescaled by a positive constant [2] . Now the horizon fields 
are 'dynamical', i.e., r-dependent, and the rescaling free- 
dom is by a positive function of r. 

We are interested in calculating the flux of energy 
associated with £^ for any radial coordinate r. De- 
note the flux of matter energy across AH by J 7 ^ '■= 
J AH T ab T a ^ b r - ) d 3 V. By taking the appropriate combina- 
tion of (1) and (2) we obtain 

T ™ =ii^J AH Nr {Hs + 2?aH ^ } d3y ■ (9) 

Since H is foliated by two-spheres, we can perform a 2 + 1 
split of the various quantities on H. Using the Gauss 
Codazzi relation we rewrite 1Z in terms of quantities on 
S: 

K = K + K 2 - K ab K ab + 2D a a a (10) 

where a a = r b D b r a - r a D b r b . Next, the fact that the 
expansion 8^ of £ a vanishes leads to the relation 

K + K = K ab r a r b . (11) 

Using (10) and (11) in eqn. (9) and simplifying, we obtain 
the result 

/ N r Ud 3 V = 16nG [ T ab ? a tf r) d 3 V 

J AH JAH 

+ [ AV{|a| 2 + 2|C| 2 } d 3 V (12) 

JAH 

where \a\ 2 = <7 ab a ab with <j ab := q a c q b d W c £ d - 
\qabq cd ^c£d, the shear of f, and |C| 2 = CCa with 
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£<j ._ q ab r c \7 c £ b . Both u ab and £ a are tensors intrin- 
sic to S. To simplify the left side of this equation, note 
that the volume element d 3 V on H can be written as 
d 3 V = N~ 1 dr d 2 V where d 2 V is the area element on S. 
Using the Gauss-Bonnet theorem, the integral of N r lZ 
can then be written as 

/ N r Kd 3 V = ( dr((fnd 2 v)=8ir(r2-r 1 ). 
Jam J ri \Js J 

(13) 

Substituting this result in cqn. (12) we finally obtain 

This is the key result we were looking for. Let us now 
interpret the various terms appearing in this equation. 
The first integral on the right side of this equation is 

(r) 

the flux Tm of matter energy associated with the vector 
field Since £?s is null and f time-like, if T a b sat- 

isfies, say, the dominant energy condition, this quantity 
is guaranteed to be non-negative. It is natural to inter- 
pret the second term as the flux Tg r ^ of £" r ^ -energy in the 
gravitational radiation: 



9 ■ 16nGj AH 



iV r {| ( r| 2 + 2|C| 2 } d 3 V. (15) 



This expression shares four desirable features with the 
Bondi-Sachs energy flux at null infinity. First, it does 
not refer to any coordinates or tetrads; it refers only to 
the given dynamical horizon H and the evolution vec- 
tor field £( r y Second, the energy flux is manifestly non- 
negative. Third, all fields used in it are local; we did not 
have to perform, e.g., a radial integration to define any 
of them. Finally, the expression vanishes in the spheri- 
cally symmetric case: if the Cauchy data (q a b, Kab) and 
the foliation on H is spherically symmetric, a a b = and 
C a =0. 

To conclude this section, let us choose for our radial 
coordinate the area radius R := a/ Air. Then, 



R2 _R]_ = wh) , wh) 
2G 2G m 9 



(16) 



Thus, as promised in part I, we have obtained an explicit 
formula relating the change in the area of the horizon to 
fluxes of matter and gravitational £(fl)-energy. 



IV. MASS AND THE FIRST LAW 

Let us now combine the results of parts II and III to 
obtain the physical process version of the first law for H 
and a mass formula for an arbitrary cross-section of H. 

Denote by E*( RS > the £(#)-energy of cross-sections S of 
H. While we do not yet have the explicit expression for 



it, we can assume that, because of the influx of matter 
and gravitational energy, will change by an amount 

A E^ R ) = Tm^ +Fg R ^ as we move from one cross section 
to another. Therefore, if we define effective surface grav- 
ity ha associated with £^ as Ur := 1/2R, the infinites- 
imal form of (16) implies (k R /8wG)da = dE^ R K For a 
general choice of the radial coordinate r, (14) yields a 
generalized first law: 



8ttG 



da = dE^ 



(17) 



where the effective surface gravity k r of is given by 

■ (18) 



dR 



where 



Mr 



(r) 



This rescaling freedom in surface gravity is analogous to 
the rescaling freedom which exists for Killing horizons, 
or more generally, isolated horizons. The new feature in 
the present case is that we have the freedom to rescale 
the surface gravity (and £ a ) by a positive function of the 
radius instead of just by a constant. This is just what 
one would expect in a dynamical situation. Finally, note 
that the differentials appearing in (17) are actual varia- 
tions along the dynamical horizon due to an infinitesimal 
change in r and are not variations in phase space as in 
some of the formulations [2, 3, 7] of the first law. 

To include rotation, pick a vector field <p a on H such 
that ip a is tangent to the cross-sections of H, has closed 
orbits and has affinc length 2tt.( At this point, (p a need 
not be a Killing vector of q a b-) Consider time evolution 
vector fields t a which are of the form t a = N r i a — tttp a 
where N r is a permissible lapse associated with a radial 
coordinate r and il an arbitrary function of r. Evaluate 
the quantity J AH T ab ? a t b d 3 V using (3) and (14): 



T2 - n 
2G 



Mi nrd ' v -i ar " 2v 



- [ 2 dQ(f f d 2 v\ = [ T ab ? a t b d 3 V 

Jfl! JS J J AH 

-tL L sip ° bc « q " >d3v ■ (19) 



Again, if we denote by E 1 the t-energy associated with 
cross-sections S of H, the right side of (19) can be inter- 
preted as AE l . If we now restrict ourselves to infinites- 
imal AH, the three terms in the curly brackets combine 
to give d(ttJ) — Jdfl and we obtain 



dr — k r 

h FldJ = — — 

2G 8irG 



da + ndJ = dE l . 



(20) 



This equation is our generalization of the first law for dy- 
namical horizons. Since the differentials in this equation 
are variations along H, this can be viewed as a 'physi- 
cal process version of the first law'. Note that for each 
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allowed choice of lapse N r , angular velocity O(r) and vec- 
tor field tp a on H, we obtain a permissible time evolution 
vector field t a = N r £ a — fl(p a and a corresponding first 
law. This situation is very similar to what happens in the 
isolated horizon framework where we obtain a first law for 
each permissible time translation on the horizon. Again, 
the generalization from that time independent situation 
consists of allowing the lapse and the angular velocity to 
become r-dependent, i.e., 'dynamical'. 

For every allowed choice of (N r , 0(r), ip a ), we can in- 
tegrate eqn. (20) on H to obtain a formula for E* on any 
cross section but, in general, the result may not be ex- 
pressible just in terms of geometric quantities defined lo- 
cally on that cross-section. However, in some physically 
interesting cases, the expression is local. For example, 
In the case of spherical symmetry, it is natural to choose 
$1 = and R as the radial coordinate in which case we 
obtain E* = R/2G. This is just the irreducible (or Hawk- 
ing) mass of the cross-section. Even in this simple case, 
(19) provides a useful balance law, with clear-cut inter- 
pretation. Physically, perhaps the most interesting case 
is the one in which q a b is only axi-symmetric with ip a 
as its axial Killing vector. In this case we can naturally 
apply, at each cross-section S of H, the strategy used 
in the isolated horizon framework to select a preferred 
t a : Calculate the angular momentum J defined by the 
axial Killing field ip, choose the radial coordinate r (or 
equivalently, the lapse N r ) such that 



and choose £1 such that 



R 4 - AG 2 J 2 
2R 3 VR i + 4G2J2 



n = n o (R) := 



2GJ 

rTWTW 2 ^' 



(21) 



(22) 



This functional dependence of k r on R and J is exactly 
that of the Kerr family. (The condition on surface gravity 
can always be implemented provided the right side of (21) 
is positive, which in the kerr family corresponds to non- 
extremal horizons. The resulting r and N r are unique.) 
With this choice of N r and fl, the energy Eg is given by 
the well known Smarr formula 



E 



2( 



k a a 
8^G 



VR 4 + AG 2 J 2 
2GR 



(23) 



as a function of its angular momentum and area, each 
cross-section is assigned simply that mass which it would 
have in the Kerr family. However, there is still a balance 
equation in which the flux of gravitational energy T^°^ is 
local and positive definite (see (19)). (The gravitational 
angular momentum flux which, in general, has indeter- 
minate sign vanishes due to axi-symmetry.) Motivated 
by the isolated horizon framework, we will refer to this 



canonical E ° as the 



associated with cross-sections 



S of H and denote it simply by M. Thus, among the 
infinitely many first laws (20), there is a canonical one: 



dM 



k Q 



da + fl dJ . 



, 8ttG 

We conclude with three remarks 



(24) 



i) Note that the mass and angular momentum depend 
only on local fields on each cross section S and changes 
in these quantities over finite regions AH of H have been 
related to matter and gravitational radiation fluxes, de- 
termined by the local geometry of H. 

ii) Unlike the vector fields — N r £ a , general permis- 
sible vector fields t a is not necessarily causal. Therefore 
the matter flux J AH T a i,t a T h d 3 V need not be positive. 
Similarly, if ip a is not a Killing field of q a b, the gravita- 
tional flux need not be positive. Therefore, although the 
area a always increase along f a , E* can decrease. This 
is the analog of the Penrose process in which 'rotational 
energy' is extracted from the dynamical horizon. 

iii) While the infinitesimal version eq (20) of the first 
law is conceptually more interesting, the finite balance 
equation (19) is likely to be more directly useful in the 
analysis of astrophysical situations. In particular, the 
presence of an infinite number of these balance equations 
can provide useful checks on numerical simulations in the 
strong field regime. 
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